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Abstract 

Planetary orbits, being conic sections, may be obtained as the locus of intersection of planes and 
cones. The planes involved are familiar to anyone who has studied the classical Kepler problem. We 
focus here on the cones. 

' p ' 1 Introduction 

Derivations of the orbit equation in the Kepler problem traditionally begin with a reduction to motion in a 
plane via conservation of angular momentum. This simple, appealing, and effective step is immediate from 
the rotational symmetry of the force law. What if it were possible similarly to "read off" from the form of 
the equations of motion that solutions must be confined to the surface of certain cones? At the least, this 
would provide a pedagogically appealing demonstration of the conic section nature of orbits, and, moreover, 
* ^ . one that is faithful to the provenance of these curves in solid geometry. 

In this paper we study the cones and cylinders of Keplerian orbits. For example, we show in section 4 
. that orbits of energy ft, < lie on the surface of either of two cylinders with axes that pass through the center 



of the orbit in the direction of 1 ± y ^^A, where 1 is the angular momentum vector, m is the mass, and A 
is the Laplace-Runge-Lenz vector. 

A conic section uniquely determines the plane in which it lies, but there are infinitely many cones which 
intersect with a plane to produce a given conic section. It turns out that the locus of vertices of all such cones 
forms another conic section in a plane perpendicular to the original plane, and with an eccentricity equal 
to the reciprocal of the eccentricity of the original conic section. We provide complete and self-contained 
proofs of these rather striking facts in section 2 below. In section 3 we recall the Kepler problem and give 
rigorous proofs of some qualitative results about orbits. Section 4 is devoted to characterizing the cones and 
cylinders of planetary orbits in terms of initial conditions, and in section 5 we calculate the location of the 
cylinder axes for the earth's orbit. 

Probably nothing in this paper is new. (How could there be anything new to discover about conic sections 
or the Kepler problem?) Modern textbooks, however, invariably view conic sections as planar curves, perhaps 
to the detriment of some very pretty mathematics that only comes into focus when these curves are viewed 
in a fully 3 dimensional context. We certainly do not claim our approach is pedagogically superior to any of 
the more standard ones, but hope that it may have some interest. 

Throughout the paper we use bold face for vector quantities. Scalar quantities, and, in particular, the 
magnitude of a vector, are given in normal typeface. Thus, the magnitude of v is v. Dots over variables 
indicate derivatives with respect to time. 
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2 Ellipses in Space 



Let E be an ellipse, considered as a point subset of 3 dimensional Euclidean space. The ellipse uniquely 
determines the plane, II, that contains it. For study of the plane geometric properties of E it is convenient 
to adopt as definition of an ellipse the following: 

Definition 2.1. An ellipse is the locus of points in the plane, the sum of whose distances from a given pair 
of distinct points F± and F2 (its foci) is equal to a given positive constant. 

(It is customary to include the circle as a degenerate case of the ellipse.) 

To state the definition of ellipse as a conic section, it is necessary to clarify some terminology and notation. 
First, by the word cone we shall always mean "right circular cone." More precisely, given a line A and a 
point on it, V, the cone C with vertex V and axis A is the locus of points P such that the smaller positive 
angle from segment PV to line A has a fixed positive value. The lines PV are called generators of C. The 
maximum plane angle, 9 G (0, 7r), between two generators is a characteristic of C that does not appear to 
have a standard name. We shall take the liberty of coining the term bloom for the angle 6. 

A cone C may naturally be divided into 3 disjoint subsets: 

C = C-U{V}UC+, 

where C± are congruent connected sets. 

If A has an orientation, designate the nappes C + and C_ in the obviously consistent way; otherwise, 
choose them indifferently 

The definition of ellipse as conic section follows: 

Definition 2.2. An ellipse is the locus of intersection of a given cone C and a given plane that is not parallel 
to any generator of C, provided the plane meets only one nappe of C. 

We shall call the cones C having E = C n II the cones belonging to E. Our goal in this section is to 
characterize the set of cones belonging to a given ellipse. 

For an ellipse, represented as the intersection of a cone C and a plane II, a Dandelin sphere is a sphere 
centered on the axis A of C, which meets one nappe of C in a circle centered on A, and which meets II 
in exactly one point. A given ellipse with accompanying cone has exactly two Dandelin spheres lying on 
opposite sides of II. 

A standard use of Dandelin spheres (see, e.g., [TJ pp. 320-321]) is to derive definition 2.1 as a consequence 
of definition 2.2, yielding along the way the fact that the points of intersection of II with the Dandelin spheres 
are the foci of E. We use them in a similar way to prove 

Theorem 2.1. Let E be an ellipse lying in a plane LT with foci F± and F 2 . The locus of vertices of cones 
belonging to E is a subset of the unique plane II' through line F1F2 and perpendicular to II. 

Proof. Let C be a cone belonging to E. Let Si and S2 be the two Dandelin spheres, which meet LT in points 
Fi and F2 respectively. Let C\ and C2 be the respective centers of these spheres. Since each sphere is 
tangent to LT, it follows that the lines C\F\ and C2F2 are each perpendicular to LT. Therefore these lines are 
parallel to each other, and since they are distinct, they determine a unique plane LT'. Since both C\ and C'2 
lie in LT', it follows that the axis of C is a subset of LT'. In particular, the vertex of C lies in LT'. The fact that 
LT and LT' are perpendicular to each other is immediate from the construction of LT'. □ 
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Theorem 2.2. Let E be an ellipse as in Theorem 2.1. Let segment AP be the major axis of E, and determine 
the plane II' as in Theorem 2.1. Then the vertices of cones belonging to E lie on the hyperbola in II' having 
foci A and P and vertices at Fx and F 2 . 

1 Proof. Let C be a cone belonging to E with vertex V. Refer 

to figure 1, in which the plane of the figure corresponds to II'. 
The circle there represents the intersection of II' with one of 
the Dandelin spheres of E. Point F represents one of the foci 
of E, and points G and H represent points of contact between 
the Dandelin sphere and C. Since circle FGH is the incircle of 
triangle APV, its center O is at the intersection of the bisectors 
of the internal angles of APV. Denoting by \AH\ the length of 
segment AH, we have \~AH\ = \~AF\, \PG\ = \PF\, and \VH\ = 
\VG\. Therefore, 

\VA\-\VP\ = (\VH\ + \AH\) - (\VG\ + \PG\) 
= \AH\ - \PG\ = \AF\ - \PF\. 




Figure 1: Intersection with plane II' 



The absolute value of the last written expression reduces to 2ea, 
where e is the eccentricity of E, and a is the semi-major axis of 
E. Since the difference of distances from all such V to the points 
A and P is constant, it follows that the locus of points V lies on 
the stated hyperbola. □ 
Denote by E 1 the hyperbola of Theorem 2.2. It has foci at A and P, one vertex at F, and eccentricity 
i. The other vertex of E' is at the other focus of E. (The excircle of triangle APV, having center at the 
intersection of VO and the bisectors of the exterior angles of APV at A and P, touches AP at the other 
vertex of E 1 , and is also the intersection with II' of the other Dandelin sphere of E.) 

Next we show that every point of the hyperbola E' , except for 
its vertices, is the vertex of some cone belonging to E. Consider 
such a point V lying on the branch of E' through F. Assume 
\AF\ > | FP |, the contrary case being similar. Construct triangle 
APV in II', and its incircle H'F'G' centered at O', as in figure 
2. 

By the same reasoning as in the proof of Theorem 2.2 
have 



we 



(2.1) 



\AF'\ - \PF'\ = \VA\ - \VP\ = 2ea. 



In particular, \AF'\ > \PF'\. 
Let C be the cone with vertex V, 



axis VO', and generators 
V A and VP. The intersection of C with II is an ellipse (by 
definition 2.2), and this ellipse has the sphere with center O' 
and radius 




Figure 2: V lying on E' 

2ae. Thus E 



|0'.F'| as one of its Dandelin spheres. It also has AP 
as major axis. Let e' be the eccentricity of this ellipse. Then it follows from (2.1) that 2ae' 
= C n II, and V is the vertex of a cone belonging to E. 

The following theorem summarizes much of the foregoing: 



Theorem 2.3. Let E be an ellipse with eccentricity < e < 1 lying in plane LT and having major axis AP. 
Let LT' be the plane through AP perpendicular to LT. Then for each < 9 < n there are exactly 4 cones of 
bloom 9 belonging to E. Each cone is obtainable from the others by reflection in plane IT and in the plane 
through the center of the ellipse perpendicular to both LT and LT'. The vertices of all such cones lie on the 
hyperbola in IT' having foci at A and P and eccentricity - . 
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The axis of a cone belonging to a non-circular ellipse intersects the major axis in a point that lies strictly 
between the center of the ellipse and the focus that is nearer to the cone vertex. Referring again to figure 1, 
in which F is the focus closest to V, we have \AF\ > \FP\, hence \VA\ > \VP\. Denoting by J the point of 
intersection of VO and AP, we have |AF| > \AJ\ > \PJ\. To see the left-hand inequality, note that, since 
the interior angle of APV at P exceeds the angle at A, the interior angle of triangle VJP at J is less than 
a right angle. On the other hand, OF is perpendicular to AP. To see the right-hand inequality, apply the 
law of sines in triangles AJV and VJP to show that the ratio of \AJ\ to \PJ\ is equal to the ratio of \VA\ 
to \VP\. 

The intersection of a cylinder and a plane not parallel to the cylinder is also an ellipse. We shall say that 
the cylinder belongs to the ellipse. The following result is the analogue of Theorem 2.3 for cylinders. It may 
be proved directly or obtained as a limiting case of Theorem 2.3 as the bloom tends to zero. 

Theorem 2.4. Let E be an ellipse with eccentricity < e < 1. Then there are exactly two cylinders 
belonging to E. The axes of these cylinders pass through the center of the ellipse and form the asymptotes of 
the hyperbola of loci of vertices of cones belonging to E. 

For completeness, we conclude by stating the analogous results for circles. If E is a circle lying in plane 
II, then there is exactly one cylinder belonging to E. Its axis is perpendicular to II through the center of E. 
For each angle 6 6 (0, tt) there are exactly two cones belonging to E. The vertices of these cones lie on the 
cylinder axis and are obtained from each other by reflection in II. 



3 The Kepler Problem 

According to Newtonian mechanics the position vector r = r(t) of a body at time t moving under the 
influence of a Coulomb or gravitational soured located at the origin satisfies the equations of motion 

kr 

(3.1) r = --3> 

where k is a measure of the force per unit mass generated by the source. For sentimental reasons we shall 
call the source "the sun" and the body "the planet." Let p = mr be the momentum of the planet, where m 
is its mass. 

As is well known, 4 integrals (conserved quantities) of the motion can readily be deduced from (13.11) : The 

2 , 

3 components of the angular momentum 1 = r x p, and the energy, h = — 

We shall be concerned here only with negative energy (h < 0) orbits. Such orbits are, as it turns out, 
generally eccentric ellipses, though there are two annoying special cases: circular orbits and catastrophic 
ones that plunge straight into the sun. The latter may be avoided by the assumption 1^0. 

Theorem 3.1. 7/1^0 and h < then we have the following a priori bounds on r and p: 

/„ I 2 ™k 

(3 - 2) 2^fc <r< W ; 



l\h\ 2m 2 k 

(3.3) 1 T<P< — i— • 

mk I 

Proof. Since 1 ^ 0, p can never vanish. Thus — = ■§— — h > —h, and the right-hand side of (3.2) follows. 

For the left-hand side of (3.2), we introduce the pivotal quantity 7 = p • r, the virial of Clausius. By the 
Pythagorean theorem, 

(3.4) r 2 p 2 = 1 2 + l 2 . 

1 The source is assumed to be rigidly fixed in space. Real two body problems can be reduced to this case. 
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Since h < we have — > #— , hence p 2 r 2 < 2m 2 kr. This, in conjunction with (3.4) yields the left-hand side 
of (3.2). 

Again, by (3.4) and (3.2) we have p 2 > ^ > m 2 k 2 ' anc ^ left-hand side of (3.3) follows. 
Finally, combining h < with the left-hand side of (3.2) we have §^ < ^ < 2 "^ fc2 , and the right-hand 
side of (3.3) follows. □ 

Theorem 3.1 and standard results on existence and uniqueness of solutions of ordinary differential equa- 
tions (see, e.g., [5J pp. 15-19]) show that, given r and po satisfying r x po / and J^- — ^ < 0, the 
system p. II) has a unique solution r satisfying r(0) = r and mr(0) = po The components of r are, moreover, 
analytic functions of t G (— oo, oo). 

If / is a function of time, we define its time average, (/) by 



(/)= lim - [ f(s)ds, 



assuming the limit exists. Also denote by T = f— and V — — ^ the kinetic and potential energy respectively. 
Our next result is a very special case of the celebrated Virial Theorem of statistical mechanics. 

Theorem 3.2. Under the same hypotheses as Theorem 3.1, both (T) and (V) exist and satisfy 

(3.5) {T ) = - l -{V)=-h. 

Proof. By Theorem 3.1 each of the functions T, V, and 7 is bounded. Thus, by the fundamental theorem of 

2 , 

calculus, (7) = 0. On the other hand, 7 = ^7 — ^ = 2T + V. Since h = T + V obviously has an average 
value, it follows that both averages (T) and (V) exist and satisfy (3.5). □ 

Remark 3.1. The virial 7 itself is the derivative of a bounded function, I — ^mr 2 , the moment of inertia. 
Thus we also have (7) = 0. 

The final result of this section will be needed in Section 4. 

Theorem 3.3. Assume 1^0 and h < 0. Then 7 vanishes for infinitely many t. 

Proof. We shall, in fact, show more: that f = for infinitely many t, which suffices, since 7 = mf • r = mrr. 
Since r is a bounded smooth function, this can only fail to happen if lim^oo r(t) = r^ exists. Now there is 
only one value for this limit that is compatible with Theorem 3.2, i.e., 

(3.6) roo = 

By energy conservation, p then also has a limit, Poo, and by (3.4) and Remark 3.1 we must have TooPoo = I- 
Together with (3.6) this implies a relation between h and /: 

™3;„2 

(3.7) h = 



2l 2 

On the other hand, by (3.4), pr > I, so h = — — > - —. A little calculus shows that the latter 

quantity has a negative minimum value of — " ;2 ° , achieved if and only if r = ~^rg ■ Thus, if (3.7) holds, r 
must be constant, and in that case r is identically zero. □ 

Remark 3.2. From the above proof, it follows that (3.7) is the circularity condition: An orbit of finite 
negative energy is circular if and only if (3.7) holds at any time t (hence, at all times t.) 
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4 Cones, Cylinders, and Kepler's Laws 



We continue to assume as in section 3 that the initial conditions for equations (3.1) are such that h < and 
1 ^ 0. If the solution of (3.1) lies for all time on the surface of a certain cone (or cylinder) C, we shall call 
C an orbit cone (resp. orbit cylinder.) 

The main goal of this section is to prove Kepler's first law - that planetary orbits are ellipses with the 
sun at one focus - by the somewhat unusual means of establishing the existence of an orbit cylinder. In so 
doing, we shall also classify all possible orbit cylinders and cones, and provide a proof of Kepler's third law. 

Let 

r 

(4.1) A = p x 1 m 2 k- 

r 

be the Laplace- Rung e-Lenz (LRL) vector. 

Theorem 4.1. The vector A is an integral of (3.1), i.e., if r satisfies (3.1), then A is constant in time. 

See, e.g., Goldstein's tex10 -3! for a proof. Also see [I] for an account of the history of this remarkable 
quantity. 

As we show next, the form and character of orbit cylinders and cones may be deduced from Theorem 
4.1, but this approach seems unfair in that it is possible to derive the complete orbit equation in a few deft 
strokes [3] from conservation of the LRL vector. Accordingly, after discussing what can be deduced from 
Theorem 4.1, we present an alternative approach that does not depend on knowing Theorem 4.1. 

Theorem 4.2. Assume ft, < 0,1 ^ 0, and that (3.7) does not hold, i.e., the orbit is not a circle. Then the 
quantity |(r + 2 A) x (1 + (3A)\ is constant in time if and only if f3 — ±-J ^TKr 

Proof. Since A and 1 are orthogonal, we have 

|(r + p 2 A) x (1 + /3A)| 2 = |(r + f3 2 A) x 1| 2 + (3 2 \r x A| 2 

= |(r + /3 2 A) x l| 2 + /3 2 |r x (p x 1)| 2 
= l 2 (r 2 + 2/3 2 r ■ A + /3 4 A 2 + /3 2 7 2 ), 



where 7 = r • p is the virial introduced in section 3 
= d£ = 2 7 (h! - ssk) and 4 

at m 1 5 at 1 V m r at 



Using ^ = ^7, = 27(2- - — ), and 4jr • A) = -^7, we find that a sufficient condition for 



r + j3 2 A) x (1 + f3A)\ to be constant is 



1 ^mk „, / p 2 mk\ 1 , „ 
0= f3 2 +P 2 \- = — + 2hf3 2 



m r \m r J m 



Conversely, the conclusion ft 2 — follows unless it happens that 7 is identically zero. But in that 
case it is easy to show that (3.7) must hold, a possibility that is explicitly ruled out in the statement of the 
theorem. □ 

For given f3, A, and 1, the locus of points r satisfying |(r + f3 2 A) x (1 + /3A)| = constant > is a cylinder 
with axis passing through the point at vector —(3 2 A offset from the sun, in the direction of the vector 1 + /3A. 
Since the orbit is confined to a plane with normal vector 1, it follows from Theorem 4.2 that non-circular 
orbits lie on the locus of intersection of a cylinder and a plane that is not parallel to the cylinder, i.e., an 
ellipse. 

Turning to an analogous result for orbit cones, let < 9 < n and set p — tan(|). Let x and y be fixed 
non-zero vectors. Then the locus of points r satisfying 



(4.2) l(r - x)Xy| = P 

V 

2 Goldstein uses k for the total gravitational force on the planet when r = 1, not the force per unit mass. Thus his k 
corresponds with our mk. 



(r - x) ■ y 

9 * 



(i 



is a cone of bloom 6 with axis passing through x in the direction of y and having vertex at x + y. 

Motivated by Theorem 4.2, we seek constants a ^ 0, (3 ^ 0, and e > so that the cone with x = qA and 
y = e(l + f3A) is an orbit cone. With these choices, equation (4.2) reduces to the form Cr 2 + Dr + E = 0, 
where 

C = l 2 (l + 2mhf3 2 ) - m 4 k 2 p 2 l3 2 
D = 2m 2 kl 2 {p 2 + a) - 2m 2 kp 2 /3{el 2 + e/3 2 A 2 + a(3A 2 - j3l 2 ) 
E = (-2a/ 2 + a 2 A 2 - p 2 l 2 )l 2 - p 2 (el 2 + e(i 2 A 2 + a(3A 2 - pi 2 ) 2 . 

(In deriving these these coefficients it is useful to note that, by (3.4), j 2 — r 2 p 2 — l 2 = r 2 ( 2m , 2fc + 2mh) — l 2 .) 
Setting C, D, and E equal to zero and solving for a, (3, and e we find that 

(4.3) e = ± ~k~ 

m^kp 



(4-5) a = - lHl 

Zmn 

In summary, we have 

Theorem 4.3. Assume h < 0,1 ^ and that (3.7) does not hold. Let < 9 < tt,p = tan(|) and f3,a and e 
be given by (4-3)-(4-5) above. Then there are 4 orbit cones: For each of the two choices of a in (4-5) there 
are two cones of bloom 6 with axis passing through qA. One has vertex at aA + |e|(l + |/3|A) and the other 
has vertex at aA. — |e|(l — |/?|A). 

In the rest of this section we attempt to derive the existence of an orbit cylinder from the equations of 
motion in the most direct and elementary manner possible. Our method is to work in cylindrical coordinates 
and seek solutions that have constant radial polar coordinate. Unfortunately the proper choice of origin and 
orientation of coordinate axes for such a coordinate system is not obvious, and will have to emerge in the 
course of the analysis. Indeed, the origin will not be at the sun, since we know from the results of section 2 
that orbit cylinder axes only pass through the sun in the case of circular orbits. 

To simplify matters, it is convenient to make some initial choices. First, we continue to assume, as we 
have above, that 1 ^ 0, h < 0, and that the orbit is not circular. Since 7 = p • r vanishes at infinitely many 
times (by Theorem 3.3), we may and do assume that we have 7 = at time t = 0. Let L be the line through 
this initial position of the planet and the sun. Let II' be the unique plane through L and normal to po, the 
initial momentum of the planet. 

Let (w, cf>, z) be the coordinates of the planet in a cylindrical coordinate system whose origin lies on L, 
whose z axis lies in plane II', and which is oriented so that the initial cylindrical angular (0) coordinate of 
the planet is zero. (For the corresponding cartesian coordinate system, II' is the x-z plane.) We use w for the 
cylindrical radial coordinate rather than the customary r since we wish to reserve r for the sun-planet vector. 
We may write the cartesian coordinate vector of the sun as S = (S x , 0, S z ). (At this point, the orientation of 
the x and z coordinate axes in II' is arbitrary. The only further requirement we shall impose for now is that 
the choices be made in such a way that the x-coordinate of ro, the initial position of the planet relative to 
the sun, and 0(0) are both positive^.) 

Let u = S x cos 4>. Because of the choices made so far, there is a relation connecting w, z, u, and r : 



(4.6) r 2 =w 2 + z 2 -2wu-2S z z + S 2 . 

3 It may seem strange to make the choice of coordinate system part of the problem, but this is not uncommon in mechanics. 
See, for example, the solution of the Kepler problem in [3] using the Hamilton- Jacobi method. 
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The equations of motion (3.1), expressed in an equivalent cylindrical coordinate form, read as follows: 



k 

(4.7) w - w<p 2 + — (w - u) = 0, 



k 

(4.8) 2w<p + wip+ -t S x sm4> = 0, 



(4.9) z + -( z - Sz ) = 0. 

We are interested in solutions of (4.7)-(4.9) having constant w = w > \S X \, since then the set w — w is 
an orbit cylinder. The gambit of locating the origin away from the sun is compensated by the simplifications 
in (4.7) and (4.8) that result from w = w = : 



(4.10) w$ 1 = ^{w-u), 



k 

(4.11) wcj) — S x sm<j). 

Indeed, equations (4.10) and (4.11) are readily integrated. Dividing (4.11) by (4.10) and rearranging, we 
have 

|ln ( 0) = -|lm>-u), 

hence 

(4.12) — ,c>0. 

w — u 

Equation (4.12), in turn, is equivalent to Kepler's Equation: 
Theorem 4.4. A C 1 function (f> satisfies (4-12) with (f>(0) = if and only if 

(4.13) <f>(t)-(^Jwi<t>(t) = lt 

holds for all — oo < t < oo. 

Substituting (4.12) in (4.10), we find that (4.10) and (4.11) entail that r must be an affine function of u: 



(4.14) r = X(w-u), A = 



k \ 3 

^27 



Conversely, if (j> = <p(t) satisfies (4.13) with 0(0) = 0, and we put r — r(t) = \{w — S x cos </>(£)), then it is 
easy to check that (4.10) and (4.11) hold. 

We have yet to consider the z equation, (4.9). Towards this end, note that (4.10) and (4.11) imply 

(4.15) u=4(— -«)■ 

r A w 

Next, by completing the square in (4.6) one has the following purely algebraic result. 
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S 2 1 

Theorem 4.5. If (4-14) an d (4-6) hold, then z is affine in u if and only if X = (1 — ;^f)~ 5 . If so, then 

(4.16) z-S z =±A(^-u). 

Equation (4.9) is an immediate consequence of (4.15) and (4.16). 

As is well known, initial values ro and po for the position and momentum of the planet relative to the 
sun determine a unique solution of the equations of motion. Thus, assuming the parameters w, S x , and the 
constant of integration c in (4.12) can be chosen in a way that is consistent with the initial conditions and 
the requirement on A from Theorem 4.5, the solution must be exactly that one given by (4. 13), (4. 14), and 
(4.16). 

Four of the degrees of freedom in the vector values of r and p have already been used to fix the line 
L and plane II'. As we show next, the remaining two degrees of freedom, determined by the magnitudes po 
and ro, together with the algebraic requirement on A in Theorem 4.5, do indeed determine one and only one 
choice of the parameters. 

There are two cases depending on whether I 2 > m 2 kr or I 2 < m 2 fcr , the case of equality being ruled 
out by the assumption that the orbit is not a circle. It turns out that the planet is initially at perihelion in 
the former case, and initially at aphelion in the latter case. We consider only the first case in detail, since 
the two are very similar. 

The constraints on w, S x , and c imposed by the initial conditions are summarized in the following three 
equations: 



(4.17) 




(4.18) I = mk3(w 2 c)3, 

and 

Equation (4.17) is a consequence of (4.14). To obtain (4.18), use / = ropo = mroW(j>(0) , (4.12), and (4.17). 
To solve these equations, it is convenient to let e = ■^ L , since this will turn out to be the eccentricity of 
the orbit. First multiply equations (4.18) and (4.19), obtaining 

I 



(4.20) (q : 

After substituting this result in (4.17), we find 

72 

(4.21) l + e = 



m 2 kr ' 

which determines the value of the quantity e from the initial conditions. (We have that e > 0, hence S x > 0, 
since the right hand side of (4.21) is greater than one - see the paragraph before (4.17). It then follows from 
(3.2) that e satisfies < e < 1.) Equation (4.19) can be solved for c using (4.20), and finally w can be found 
from (4.20), (4.22), and (4.17.). The results are: 

(4.22) c= mfc(1 - e2) , 
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(4.23) w= o^)m, 

1 — e 

and 

(4.24) S x = we = er 



1 + e 

1 - e 



(In case 2, one takes e = — ^ and proceeds similarly.) 

It follows from similar triangles that ^- — r "g S , from which we conclude that 

(4.25) S = r 



1 -e 



The location of the origin is now completely determined. It is on the line L at distance S from the sun in 
the direction opposite the planet. To complete the description of the coordinate system, choose the positive 
z direction so that the coordinate system is right-handed. Then the z-component of 1 is positive, which 
implies that the angle, ip, from the sunward side of L to the positive z-axis satisfies < ip < ir. Since one 
has simp = ro +g i it is easy to derive from (4.23) and (4.25) the following rather pleasing result: 

(4.26) e=\cos(ip)\. 

The two solutions for ip correspond to two orbit cylinders. If ip is acute then both S z and the z-component of 
r are positive. If ip is obtuse, then both are negative. In either case, another application of similar triangles 
yields 



(4.27) | 5z | =ro _^_. 

Letting a = j^, we may write (4.14), the orbit equation, as 

(4.28) r = a(l - ecost^), 

where <p> is determined by solving Kepler's equation (4.13). This form reveals that the angular cylindrical 
coordinate <p may be interpreted as the eccentric anomaly, and that the radial cylindrical coordinate, w, is 
equal to b, the semi-minor axis of the orbit. 

If < -0 < | then S z > 0. If we take the minus sign solution from (4.16) we obtain 

(4.29) z = aecoscj). 

(We can eliminate the positive sign solution because, by (4.25) and similar triangles, we have z(0) = ae.) If 
^ < ip < n then S z < 0. In this case, taking the positive sign solution from (4.16) yields z = —aecos(p. 

It is worth noting the character of the motion in various 2 dimensional projections other than the orbital 
plane. The pair (x, y) moves around the circular cross section of the orbit cylinder: (x, y) = (w cos <j>, w sin <p). 
The motion of the pair (x, z) is degenerate, being confined to the line L through the origin of slope —j==^, 

while the pair (z,y) moves on an ellipse centered at the origin with axes of lengths ea and y/l — e 2 a. 

We conclude this section with a derivation of Kepler's 3rd law. Let t be the period of the orbit. Then, 
by (4.13) we have r = 2?rf . Thus, by (4.20) and (4.21), 

2irla 



rri 



kVT 



To eliminate I, use Kepler's 2nd law (an easy consequence of conservation of momentum) to find It = 2nabm, 
thus obtaining Kepler's 3rd law: 



2 471-2 3 



10 



5 The Orbit of the Earth 

The orbit cylinder axes, or equivalently, the asymptotes of the set of orbit cone vertices, may be regarded 
as elements of an orbit. It is therefore of interest to locate these axes for significant orbits such as the orbit 
of the earth. 

The orbit of the earth is not a fixed ellipse, so it is necessary to begin by clarifying what we mean by 
"orbit of the earth." We shall regard the "earth" as an idealized point representing the center of mass of the 
earth-moon system. Then, following the usual procedure for reducing a two body problem to a one body 
problem, we regard this point as being in orbit around the center of mass of the solar system, rather than 
the center of mass of the sun. The resulting orbit is still not an ellipse due to the influence of the other 
planets. Nevertheless, during a reasonably short time interval - a year, say - the orbit is an ellipse to very 
high precision. Accordingly, we shall define the elliptical orbit of the earth as the best approximating ellipse 
during a selected year. It is natural to select for this role the year 2000, since modern tables of stellar and 
planetary positions are referenced to the ecliptic and vernal equinox as they existed at 12 hours coordinated 
universal time on Jan 1, 2000 (Julian date 2451545.0, or J2000 for short.) The Jet Propulsion Laboratory web 
site [5] provides the following values for selected elements of this orbit: semi-major axis a = 1.00000261AU, 
eccentricity e = 0.01671123, inclination of orbital plane to reference plane I — — 0°. 00001531, and longitude 
of the perihelion it = 102°. 93768193. (Recall that the longitude is measured eastward along the ecliptic from 
the J2000 vernal equinox.) 

The nonzero value of the inclination element deserves a comment. Tabulated celestial coordinates refer- 
ence the mean ecliptic and equator at the epoch J2000, not the actual ecliptic and equator. The mean in 
question is an average over periodic perturbations ( e.g., nutation of the earth's axis of rotation,) which are 
more significant for the equator and equinox than for the ecliptic. Since we are, in any case, considering an 
idealized orbit, we shall simply assume the orbit lies in the mean ecliptic and take / = 0.0. 

The orbit cylinder axes intersect the celestial sphere at 4 points which may be regarded as lying infinitely 
far from the observer. Thus, the geometric and apparent positions of these points coincide. (Precise cal- 
culation of apparent positions of real celestial bodies - planets, moons, and even stars - must address the 
location of the observer, effects due to motion of the observer, and the finite speed of light.) We shall identify 
the axes by their points of intersection with the celestial sphere. 

Since the cylinder axes lie in a plane perpendicular to the ecliptic through the perihelion, the circle of 
celestial longitude (A) of all 4 points is the same as for the perihelion. The celestial latitude, j3, is the same as 
the angle we denoted ip in (4.26) above, and is determined by the eccentricity alone. The celestial coordinates 
of the 4 points are then as follows, where we list the pairs (/3, A), and in brackets the corresponding declination 
and right ascension: (89°. 0425, 102°.9377)[+67°.4930, 17h57m45.6s], 

(89°. 0475, 282°.9377)[+65°.6270, 18h2m4.7s],(-89°.0425, 282°.9377)[-67°.4930, 5h57m45.6s], 
and (-89°. 0425, 102°.9377)[-65°.6270, 6h2m4.7s]. 
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